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SYMBOLS ‘ 


c 

Cp 

s 

d, D, Oj 
L 
iL 
M 

^.25 

R 

f^D/Sp 
X, y, z 


configuration chord length 

pressure coefficient 

Interference pressure coefficient 

jet diameter 

aerodynamic lift 

interference aerodynamic lift 

aerodynamic pitching moment 

interference pitching moment about the jet center 

interference pitching moment about the configuration 
quarter chord 

jet-to-crossflow velocity ratio 

ratio of the jet diameter to the control point 
spacing 

axes of the jet coordinate system^ see Figure 1 


K, Y, Z axes of the panel code coordinate system, see Figure 2 

XdCT, YJCT» ZilCT location of the jet center In the panel code 
coordinate center 

Uj jet exit velocity 

freestream velocity 

a model angle of attack 

S model yaw angle 

Y angle between the jet z>ax1 s and the panel code 

Z-axIs, jet roll angle 


ill 
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jet Injection angle 
Isentropic jet thrust 



summarV 


A method Is proposed to combine a numerical description of a let In 
a crossflow with a lifting surface panel code to calculate the 
jet/aerodynamlc-surface Interference effects on a V/STOL aircraft. An 
Iterative technique Is suggested that starts with a model for the 
properties of a jet/flat plate configuration and modifies these 
properties based on the flow field calculated for th^s configuration of 
Interest. The method would estimate the pressures* forces, and momeitts 
on an aircraft out of ground effect. 

A first-order approximation to the method suggested Is developed 
and applied to two simple configurations. The first-order approximation 
is a non-iterative procedure which does not allow for Interactions 
between multiple jets In a crossflow and also does not account for the 
Influence of lifting surfaces on the jet properties. The jet/flat plate 
model utilized In the examples presented Is restricted to a uniform 
round jet Injected perpendicularly into a uniform crossflow for a range 
of jet-to-crossflow velocity ratios from three to ten. Numerical 
results for a streamlined body of revolution and a symmetrical airfoil 
are presented. The numerical fesult$ show that there Is good agreement 
between experimental and model calculated surface pressure data for the 
body of revolution and the non-lifting wing, but Indicate the need for 
iterative techniques for handling the Interactions between lifting 
surfaces and a jet In a crossflow. 

V 



This report Is divided Into two volumes* The first volume Isa 
theoretical description of the computer code developed* The second 
volume Is a detailed Users Guide for the computer code* 




SECTION I 
INTRODUCTION 

Vertical or short take-off and landing (V/STOL) aircraft 
traditionally supplement aerodynamic lift Mith high energy lift Jets* 
While the use of lift Jets Is efficient when viewed from a static 
standpoint* they create adverse aerodynamic Interference effects when 
used to supplement wing-generated lift during forward flight. When lift 
Jets are used during forward flight they are called jets in a crossflow, 
and the mode of flight is termed transitional flight. The adverse 
interference effects result from the viscous mixing of the lift jet 
fluid with the crossflow fluid. The major adverse effects are a lift 
loss and a nose-up pitching moment. These effects increase in severity 
with increasing crossflow velocity. 

As defined herein a jet in crossflow refers to the qualitative 
features of a jet injected into a crossflow. The terminology 
jet/aerodynamic- surface refers to a jet in crossflow issuing from an 
arbitrary aerodynamic surface. A jet/fl at plate refers to a jet in 
crossflow issuing from a flat plate configuration. 

Modern high speed computers, through the use of panel codes, can be 
used to solve for the flow field on or about a conventional aircraft for 
the conditions of attached flow. A panel code is a program which solves 
a linear partial differential equation numerically, by approximating the 
configuration surface with a set of panels, on which unknown singularity 
strengths are defined. The application of boundary conditions at a 
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discrete set of points, such as the paneV centroids, generates a system 
of linear algebraic equations relating the unknown singularity strengths 
to the flow at the panel centroids. The equations are then solved for 
the singularity strengths, which once known, provide the properties of 
the flow field about the configuration.^**^ 

^'vAr'rgh existing panel codes can be used to solve the flow field 
about a convent*! onal aircraft, they fall to predict the adverse 
Interference effects on a V/STOL aircraft in transitional flight. One 
method for predicting the |et*1n-crossf1ow Interference effects on a 
V/STOL aircraft Is that of Kuhn,^ The method developed by Kuhn applies 
empirical corrections to the aerodynamic coefficients for aircraft 
without jets In crossflow. Kuhn's method estimates the Integrated 
Interference effects; lift loss, and pitching moment. There are methods 
currently under development which will combine programs for predicting 
jet/flat plate Interference effects with panel codes.® 

Most jet-in-crossflow models In use today predict the flow field 
due to the properties of a jet/ flat plate model. This Is because the 
develo(xnent of a general jet/aerodvnamic- surface model Is too complex. 
Two programs In use today, which estimate jet/flat plate interference 
effects, use potential flow singularities to model the jet properties. 
The first, developed by Mooler et al.,® combines a distribution of sinks 
and doublets along a calculated jet path to obtain the velocity field 
due to the jet interference. The second program was originated by 
Oietz^ in his masters thesis and later developed to its present form by 
Fearn,® The Oletz/Fearn model assumes two contrarotating vortices to be 
the dominant flow feature, thus the vortex properties are the major 
Influence In the jet interference effects. A non-potential flow program 
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has been developed by Adler and Baron,^ which uses an integral control 
volunie method to represent the axial flow Internal to the jet plume* 

A jet/aerodynamfc-surface modeling method Is proposed in this 
report which combines a jet/flat plate model with a low-order panel code 
to model the flow field about a V/STOL aircraft when operating out of 
ground effect. The method proposed should be capable of accounting for 
all interactions between jets in crossflow and a lifting surface 
configuration. A jet/aerodynamic-surface model is presented! based on 
the proposed method, which accounts for the influence of jets in 
crossflow on the lifting surface, but does not account for the Influence 
of the lifting surface on the jet properties. The model presented 
combines the jet/fl at plate model of Fearn® with a low-order panel code 
developed by Analytical Methods Incorporated.^^ 

.Section U of this report discusses modeling methods, including 
panel codes and jet-in-crossflow codes. Sections III and IV deal with 
the mechanics of combining the two computer codes used and the testing 
of the resultant computer code. Section V presents the results for the 
configuration studies and concludes Volume 1 of this report. Volume 2 
of this report presents a detailed Users Guide for the UBWJAS computer 
code developed in Volume 1. 

Financial support for this study was provided by NASA Grant NCC2-133 
under the technical direction of Mr. D.G. Koenig, NASA Ames Research 
Center. This report is the result of Mr. K.L. Furlong's masters thesis, 
supervised by Or. R.L. Fearn. 



SECTION M 
MOOELIKQ METHODS 

Low-Order Panel Codes 

A panel code Is a program which solves numerically the linear 
partial differential equation describing potential flow about a 
surface* It does this by approximating the surface with a set of panels 
on which singularities of unknown strength are defined. The common 
singularities are sources to describe a non-lifting surface, and In 
addition, doublets or vortices, to describe a lifting surface, A Vow- 
order panel code is a program In which the singularity strength 
distribution is piecewise constant over the configuration. Higher-order 
panel codes require the singularity strength distribution to fit linear, 
quadratic or higher-order equations* 

A generic wing-body configuration is shown in Figure 1, which is a 
reproducti on from reference U. The figure shows source panel s defi ni ng 
the configuration shape; the panel control points are shown as dots at 
the panel centroids* The vortex system is shown in the wings, with 
trailing vortices to represent the wake region downstream of the 
wings. Each •singularity shown has associated with it a perturbation 
velocity flow field. The vector sum of. all perturbation velocities and 
the freestream velocity describe the flow field of a particular 
configuration* Analytical expressions for the perturbation velocity 
fl ow fi eld i nduced by a constant Source di stri buti on On an arbi t rary 
panel are given by Hess and Smlth.^ Similarly, the velocity field 
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figure 1. - Source and Vortex Panel Arrangements on a 
Uing-Body Configuration (reference 11). 
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1ndud(^%’the,«le(netits of a vortex distribution are given by Woodward^ 

■ Y*f’ 3"' •-’ * 

and by Rubbert et a1 

The perturbation velocities are used to calculate the coefficients 
of a system of linear algebraic equations, relating the magnitude of the 
velocities at the panel control points, to the unknown singularity 
strengths. The singularity strengths, which produce the perturbation 


velocities necessary to satl sfy the boundary conditions on each panel, 
are determined by solving the system of linear algebraic equations. 
With the singularity strengths known, the velocity field Is determined 
and the pressure coefficients calculated using Bernoulli's equation. 
Ths pressure distribution Is Integrated numerically over the 
configuration to yield the forces and moments on the configuration. 


The panel inet hod chosen for this project Is a low-order panel code 
composed by Dvorak, Ma skew, and Woodward.^^ The method Is known by the 
acronym WBABL, which stands for Wing Body Aerodynamic program with 
Boundary Layer solutions. The major reasons for choosing the WBABL 
panel code are Its accessibility and the availability of knowledgeable 
program users at the NASA Ames Research Center, where this project 


originated. 

Jet-1 n-Crossfl ow Codes 

A jet/aerodynamic-surface model should be capable of predicting the 
Interference effects on an aerodynamic configuration due to the flow 
field associated with a jet in a crossflow. As a prelude to discussing 
jet/aerodynamic-surface models, a qualitative description of the flow 
field Is In order. 
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Oe^Hotion of the FI w Field 

Based on extensive expeHcnental investigations of a subsonic jet 
Injected Into a subsonic crossflow, a qualitative description of the 
flow field can be made. This description of the flow field follows that 
presented by Fearn.^ The characteristics of a jet In a crossflow 
described here apply mainly to a round jet. Injected perpendicularly 
Into the crossflow, from a flat plate* Although the characteristics 
given below are for a particular jet/aerodynamlc-surface configuration, 
they are characteristic of the general jet'-in-crossflow flow field. 

Jet c,c«ren fur a relatively short distance from the jet exit plane, 
the core flow of the jet Is characterized by slow changes In flow 
properties such as velocity profile and turbulent Intensity, For a 
submerged jet, with no crossflow, the jet core Is conical In shape and 
extends approximately six jet diameters from the jet exit plane^^ before 
It Is eroded Into a highly turbulent flow* When the jet Is Injected 
Into a crossflow, the length of the jet core Is found to decrease with 
decreasing jet-to-crossflow velocity ratio* 

Shear layer* The boundary between the jet and the crossfl ow can be 
considered a shear layer between the high and low energy flows* This 
shear layer Is thin near the jet exit, but rapidly thickens as the jet 
core diffuses* It Is thought that this shear layer can be represented 
by a region of concentrated vortlclty,^^ the diffusion of which cart be 
used to describe the flow field of the jet core with no crossflow* To 
the presence of a crossflow the shear layer (vortex Structure) Is 
distorted as well as diffused by the crossflow* The distortion of the 
shear layer is thought to be the origin of the pair of contrarotating 
vortices which are characteristic of the jet-in-crossflow fl ow flel d. 
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Contrarotating vofttx pair , the dominant feature In the jet-1 n- 
crossflow flow field Is a pair of diffuse contrarotating vortices* The 
vorteK pair has been observed to start near the jet exit and extend far 
downstrean*^^ The diffuse vortex pair Is deflected downstream along 
curved paths. These paths lie on either side of the plane of flow 
syninetry and are located on the concave side of the jet 
centerllne.^^*^^ As the contrarotating vortices are swept downstream 
they diffuse rapidly until their core radius Is approximately equivalent 
to the half spacing between the vortex centers. The diffusion of 
vorticity between the two vortices results In a decrease In the strength 
of each vortex. As an Indication of the dominance of the contrarotating 
vortices as a flow feature, the vortex pair has been observed as far as 
1000 jet diameters downstream from the jet orifice.^® 

Jet Plume. The jet plume is readily observed by flow visualization 
techniques such as smoke Injection. The deflection and decay of the 
initial jet can be detected by velocity measurements In the jet plume. 
The curve that traces the locations of the maximum jet speed, from the 
terminus of jet core through cross sections of the jet plume* Is called 
the jet centerline. The location of the jet centerline can be 
determined to a point where the local maximum of the jet velocity 1s 
Indistinguishable from the free stpeam velocity. The extent of the jet 
centerline has been found to be fifteen to twenty jet diameters 
downstream of the jet exit for the Instrumentation used in the tests of 
Fearn and Weston^^ (1.e.» yaw/pitch probes). 

Wake, implicating the flow field near the jet exit Is a wake 
region. This wake Is clearly visible In oil smear studies on a flat 
pi ate. The wake region begins just downstream of the jet core and is 
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a result of a separation of the crossflow as It flows around the jet 
core* 

Most of the features described above are shown In Figure 2* The 
figure Is representative of a jet In a crossflow Issuing from a flat 
plate at a ninety degree Injection angle and a jet-to-crossflow velocity 
ratio (R) of eight* The stippled area represents the observed smoke 
plume, within which the relative locations of the jet centerline and 
vortex curve are shown* The vortex curve Is a projection of one of the 
vortex paths Onto the symmetry plane* The diffuse contrarotating vortex 
pair Is Illustrated In a cross section to the jet plume, and relative 
positions of the jet core and flat plate wake region are shown. Also 
shown In Figure 2 are the coordinate axes used to describe the jet 
properties. 

Jet/Flat Plate Models 

A jet/aerodynamic-surface model should provide the perturbation 
velocities at a point of Interest due to the combined Influences of the 
jet/aerodynamic-surface structure* Existing modeling techniques for 
transitional flight consist largely of empirically derived models, using 
potential flow singularities to account for the complex flow 
structure* Three jet/aerodyhamlc^surface models are discussed here, all 
use a flat plate as the aerodynamic-surface configuration. 

The first model discussed, developed by VIooler et a1*,^ combines 
distributions of sinks and doublets along the calculated jet path to 
obtain the velocity field due to the jet Interference* As shown In 
Figure 3, the sinks are uniformly distributed along an axis normal to 
the freestream, and account for the entrainment of the jet. Their 
strength, which varies with the distance from the jet. Is dependent upon 




Figure 2. - Prominent Features of a Jet in a Crossflow. 
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three MpirlcaVIy determined coefficients. The values chosen result In 

# 

a good correlation between experimentally and theoretically derived jet 
centerline locations. The doublets are distributed along the jet 
centerVIne and account for the jet blockage effect by creating a flow 
past an equivalent circular cylinder. This model Is effective In 
predicting the Induced perturbation velocities for aircraft 
configurations which do not include Surface area behind a jet In 
crossflow. For configurations which do Include surface area behind a 
jet In crossflow, Uooler's model, along with others, underestimates the 
magnitude of the negative pressures found in the wake region. 

The Wooler model has been extended to 1 nclude a wake correction by 
Walters and Yen.^^ This extension to the Wooler model consisted of a 
pair of Short-lived, low-strength vortices Incorporated Into the wake 
area. The addition of these vortices resulted in an improvement In the 
prediction of the surface pressures on the flat plate. While this wake 
model improved the ability of the Wooler model to predict the pressure, 
force, and moments bn the flat plate model, it still was not able to 
match the negative pressures found Immediately aft of the jet core. 

This wake model, though applied to the Wooler model, could be adapted to 
other jet/aerodynamic-surface models. 

An extension of the Wooler model has been made to model multiple 
jets In crossflow. This work, presented by Ziegler and Wooler, gives 
methods for modeling both tandem and transverse jet/flat plate 
configurations. The model does not account for any viscous wake 
effects. In thi s method the case of transverse jets 1 s modeled by 
assuming independent jets until the trajectories meet, then combining 
the jet plumes. For the tandem jet case, the crossflow speed on the aft 
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jet 1$ decreased as a result of forward jet blockage . and the two Jet 
pi urnes are combined as they meet* 

The second model was presented by Oletz^ In his masters thesis and 
later developed to Us present form by Fearn*^ This model assumes the 
two contrarotating vortices to be the dominant characteristics of the 
Jet/flat plate configuration* A less dominant* yet still Important 
characteristic* Is the entrainment of freestream fluid Into the jet 
plume. The model Is Implemented by placing a series of linear vortex 
filaments along the vortex paths* and a series of linear source segments 
along the jet centerline curve. The strength distribution for the 
vortex filadients Is determined by inferring the vortex properties from 
the experimental data of reference 14* 

Similar to the Mooler model, the model of Fearn Is unable to 
predict the pressure distribution In the wake region* Fearn *s model 
does, however, contain an empirical wake correction* which when used 
results in good correlation of model and experimental data for the 
pressure, force* and moment data for flat plate models. Fearn* s model 
has been applied to the cases of a round jet/flat plate configuration 
with perpendicular jet Injection and jet-to-crossflow velocity ratios 
from three to ten. 

The third and final model discussed In this report was developed by 
Adler and Baron.® Their model Is based on a control volume method. 

This model Is a quaslUhree-dlmenslonal Integral method used to solve 
the problem of an Isothermal, incompressible, turbulent, rourid jet In a 
crossflow. The model Is based on two Integral momentum equations, one 
written for a direction parallel to the jet centepline and the other for 
a direction perpendicular to the centerline. The mathematIcaT model Is 
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numerically solved* producing a flovi field within the jet plume* 

Results of this method, as reported in reference 9* show that It Is 
capable of predicting the axial flow In close proximity to the jet 
plume, but does not adequately predict the velocity and pressure 
perturbations on cats'lguratlons far from the jet plume* The Adler and 
Baron method Is not restricted to a jet/flat plate configuration* 

The choice of jet/ flat plate «iwdel for this project Is the method 
of Fearn. This Jet/ Aerodynamic-Surface Interference (JASl) program Is 
known as JAS IB In Its present form and Is documented in reference 8* 

Combined Methods 

There is a need for a computer modeling method which combines a 
panel code with a jet/aerodynamic-surface model* At this time one 
method Is In use which combines jet/flat plate models with panel 
codes* It is the one developed by Beatty and Kress* At the time of 
this writing there are a number of studies ongoing which are addressing 
this problem*^ Unfortunately, reports on these studies are not 
available for description here. 

The computer code developed by Beatty and Kress^- Is Intended to 
model virtually all aspects of V/STOl flight. Their computer program 
uses the Hess potential flow code to model freestream effects on the 
ai rcrafti, Suppl ementlng the Hess code 1 s a series of programs to model 
the propulsive effects on the aircraft* The propulsive effects section 
Include; an Inlet analysis secti bn, three jet/flat plate sections, and 
a ground effects section. The three jet/flat plate sections are a 
modi fled Wool er^ method, modi fled Meston^^ method, and the Thames and 
Weston^^ method for rectangular jets In crossflow. The modified Wool er 
method Is essentially the same as that described earlier except for the 



addition of a ground effects algorlthni and a limitation that there be no 
more than two jets per jet system* The Weston^^ method Is a version of 
the early work of Dietz and Pearn* The modified Weston model of Deatty 
and Kress Is reported In reference 21 to be applicable to various 
Injection angles* various jet-to-crossflow velocity ratios, and to 
multiple jets* The extension to handle multiple jets Is similar to that 
presented by Ziegler and Wooler, as mentioned earlier* The ability to 
model various Injection angles and velocity ratios Is accomplished by 
extending the work presented by Oletz^ In his masters thesis* The work 
of Thames and Weston^^ on a rectangular jet Is limited to a jet with an 
aspect ratio of four* The jet Is oriented with the major axis either 
parallel or perpendicular to the crossflow* The method used by Thames 
parall els that used by Dietz for the round jet case* 

Results for the programs of Beatty and Kress, as given in 
reference 21, show that the modified Weston model gives surface pressure 
results which agree closely with the Wooler method* Both the modified 
Weston and Wooler methods give good agreement with experimental surface 
pressure data for regions exterior to the jet wake area; the need for a 
wake model for the jet wake region 1s stressed. Another area of concern 
presented In reference 21 Is that of viscous and non-viscous Interaction 
between the various elements of the entire program package, i.e*, the 
Interactions between the freestream flow section and the propulsive 
effects section* 

Proposed Method 

The objective of this proposed method Is to determine the inviscid, 
potential flow field on and about a V/STOL aircraft In transitional 
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flight out of ground effect. The proposed method Is based on the 
following assumptions: 

1. A potential flow model, for the flow field Induced by the 
jet plume of a known jet/aerodynam1c*5urface 
configuration exists. This model Is referred to as the 
jet/base configuration model* The jet plume properties 
are assumed to be dependent on the following: the 
aerodynamic-surface configuration, jet orifice shape and 
velocity profile; jet Injection angle, and the jet-to- 
crossflow velocity ratio. 

2. The V/STOL aircraft configuration of Interest Is modeled 
with a potential flow panel code. 

3. A local velocity ratio for the jet plume properties can 
be defined. Fora jet/flat plate configuration. It is 
conventional to define the velocity ratio as either 

or Uj/U.* Where Uj is the jet exit velocity 
magnitude, and U^ is the crossflow velocity magnitude. 
For a local velocity ratio, the crossflow velocity Is 
replaced with the local onset velocity, which can be a 
vector or a scalar quantity depending on the jet plume 
model. The local onset velocity 1$ considered to be the 
velocity (or speed) at a point of interest in the jet 
plume (i.e., a singularity control point) due to all 
influences other than those of the Jet plume under 
consideration. 
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The proposed method Is Implemented by following these steps: 

1. Calculate the panel code singularity. strength 
distribution. The onset flow to the panel control points 
Is the sum of the freestream flow and all jet/base- 
configuration Induced flows. 

2. Use the results of step 1 to calculate the local velocity 
ratio for each point of Interest In the jet plume(s). 

3. Use the local velocity ratio for each point of Interest 
In a jet plume to modify the jet plume properties. 

4. Check for model convergence by comparing the jet plume 
property changes with a preset tolerance. 

5. If the model has not converged then return to step 1 with 
the new jet plume properties. If the model has converged 
then the potential flow solution Is complete. 

The above method should, after convergence, produce a model which 
represents the jet/aerodynamic-surface configuration of interest better 
than If the method were stopped after step l. The only criterion 
mentioned above for adjusting the jet plume properties Is the local 
velocity ratio. There may be other criteria which are Important to the 
jet plume properties. 

Viscous flow modelling Is not Included In the proposed method. If a 
viscous flow model. 1.e.» a wake model. Is required It could be added 
after the potential flow calculation section. 

One should note that the above method describes the modeling of 
multiple jets In crossflow. The flow field of multiple jets In 
crossflow, all in close proximity to each other, has been observed, by 
Braden, 2^ to be more compl ex than might be expected. Two better 
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procedures for handling multiple jets In crossflow nrtght be: 1) to use 
a jet/base configuration which models multiple jets In crossflow; 

2) construct a con^oslte model similar to that of the Ziegler and Wooler 
method 

Presented Hadel 

AUhough the method proposed above to solve for flow about V/STOL 
aircraft Is feasible at this time* It Is not attempted for this 
report* As a test and demonstration of the method proposed, a 
f1rst*order model, fabricated after the method described above. Is 
developed, tested, and applied to two specific configurations* The 
first-order model Is essentially the first step of the method proposed 
above* The model utilizes an empirical jet plume model from a jet/flat 
plate configuration, for the jet Influences on the lifting surface 
configuration* The model has no Interaction between multiple jets In 
crossflow, although It will accommodate multiple jets In crossflow on 
the lifting surface configuration. Also, due to the non-iterative 
design of the model, the lifting surface has no Influence on the jet 
properties* The first-order model Is limited by the jet/flat plate 
model to the case of subsonic uniform round jets Injected 
perpendicularly Into a subsonic unifom crossflow, with iet-to-crossflow 
velocity ratios from three to ten* 

The computer model presented here Is assembled by combining a 
modified version of the Jet/ Aerodynamic Surface Interference program 
(JASIB) of Fearn,^ as a set of subroutines, with a low-order panel code 
for three-dimensional lifting surfaces In subsonic Inviscld flow 
(WBABL).^® The model presented Is verified by modeling a flat plate^® 
with a rpund jet In a crossflow. Two other configurations modeled are a 



symmetrical airfoil^® and a streamlined body of revolutlon.^^ The 
results of each configuration tested are compared ytith experimental data 
In the form of surface pressure coefficients and Integrated body forces 
and moments* where they are available* 



SECTION II I 

COMBINING TWO COMPUTER COOES 

The two computer modeling codes to he combined are WBA6L* a low- 
order panel code* and JAS1B76, a jet/flat plate modeling code* The 
JASIB76 code Is a direct extension of the OASIS eodei developed by 
Fearn*^ The JASIB code models the flow over a flat plate with a round 
jet Injected perpendicularly Into the freestream flow. The JASIB76 code 
Is primarily designed to function as a group of subroutines on a 
CDC-7600 computer* and to be Independent of any surface configuration* 
The WBABl panel code Is a low-order panel code for solving the flow 
about a wing/body eomblnatlon* The WBABL code also has a boundary layer 
Iteration section which Is not utilized In the combined computer code* 
Both computer codes are written In FORTRAN IV. 

WBABl Computer Code 

The general structure of WBABL Is shown In Figure 4* the WBABL flow 
chart* The actual program Is made up of a series of overlays* to meet 
the memory requirements of the COC-7600 computer system. The overlay 
structure Is shown In Figure 5* The executive program* WBABL, shown as 
Overlay (0*0), controls the overall panel code by calling in turn the 
potential flow overlay, WBOLAY, and the boundary layer overlay, 

INTEGRAL. The calculation sequence for the WBABL panel code Is as 
follows: 

Loft geometry. AM Input parameters for the lifting and non- 
llftlng surfaces are read In the overlay WBPAN. The planar panels are 




t 



Figure 4. - WBABL Prqgrani Flow Chart. 
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lofted and the lifting surface vortex network Is laid out. The X, z 
coordinates of the four panel corner points* for all surface and vortex 
panels being considered* are calculated and written to an auxiliary file 
and the output file. 

Calculate potential flow . Within the WBAERO overlay all potential 
flow calculations are made. Due to the numerous computations made In 
this overlay* It Is the largest overlay of the entire panel code 
structure. WBAERO executes as follows; The corner points computed In 
WBPAN are read from an auxiliary file. Subroutine ELEHEN is then called 
for each panel. This subroutine calculates the control point locations 
and the transformation matrix* and then transforms the corner points 
from the reference coordinate system to the panel coordinate system. 

The panel control and corner points are then transformed to an analog 
body by subroutine ANALOG. This is required by Gothert's 
compressibility rule and Is done whenever the freestream Mach number Is 
greater than one-tenth. Subroutines INFLU and WINGIN calculate the 
velocity Influence coefficients Induced by unit source singularities and 
unit vortex singularities* respectively. These 1 nfluence coefficients 
are combined to form the aerodynamic Influence coefficient matrix. The 
boundary conditions* i.e.* zero normal flow at each panel control point, 
are calculated for each angle of attack and yaw. With the above 
information the system of algdbraic equations relating the singularity 
strengths and the boundary conditions Is solved for the source and 
vortex strengths in subroutine SOLVE. With the singularity strengths 
known, the velocities at the control points are determined. With the 
yelocities known, the surface pressures are determined. The forces and 
moments on the body are found In subroutine FORMOM, A more detailed 
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description of the above procedure Is given In the theory section of 
reference ll. 

Streamline calculation . Overlay STREAM calculates the trajectory 
of up to twenty- five body streamlines* The starting location of the 
streamlines can be any control point* The streamline Is traced fore and 
aft to the body limits. 

Boundary layer calculation. Overlay 1NTE6RAI contains the boundary 
layer calculation subroutines* which determine laminar and turbulent 
boundary layer development along the streamlines specified above* A 
detailed description of the boundary layer subroutines is given In 
reference 10* They are not utilized for this study and are not 
discussed here. 

Modification of the UBABL panel code for inclusion of a jet-in- 
crossflow module is necessary If the effects of jets In crossflow are to 
be evaluated* The Ideal solution to the problem of how to Integrate the 
jet-1 n-crossflow module into VIBABL would be to create a new overlay for 
the jet effects. This Is not possible* however, because of the overlay 
hierarchy and the eKi sting structure of VIBABL* The input to the jet 
code requires that it be addressed after the control points are defined* 
e*g.* within UBAERO* The output from the jet code and the jet 
perturbation velocities* must be available to WBAERO for 1 nclusion In 
the control point boundary condition calculation* These requirements 
dictate that a jet overlay be addressed from within the WBAERO overlay, 
but a second level overlay* (I.e., WBAERO) can not call other overlays* 

The fact that the jet-in-crossflow module can not be an overlay 
necessitates that it become a group of subroutines within the WBAERO 
overlay. This requirement brings about new complications , in the form 
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of available memory. MBAERO is tbe largest overlay of the entire WBABl. 
package, thus occupying the maximum memory. The maximum number of 
panels (control points) that a configuration can have is limited by 
ayailabU memory, which Is controlled by the size of the programs and 
subroutines within the memory. Thus, in order to have a jet-in- 
crossflow module within the panel code, the panel code must sacrifice 
its most valuable asset; its ability to model large complex 
configurations and still maintain accuracy with a fine panel network. 

JASI576 Computer Code 

The JASIB76 program is a direct descendent of the JASIB program 
developed by Fearn,® The JAStS program calculates the pressure 
distribution induced on an infinite flat plate at zero angle of attack, 
due to the influence of a single round jet in a crossflow. The flat 
plate is modeled by the method of Images, and the pressure coefficients 
at the control points on the plate are found by the application of the 
Incompressible Bernoulli's equation. The JASIB76 program is a group of 
subroutines, controlled by an executive subroutine, which computes the 
perturbation velocities, at a specified point, due to the Jet plume 
properties of a jet/flat plate model. The properties of the jet/flat 
plate model used in JASIB76 are those of a round, uniform, subsonic jet 
injected perpendicularly Into a uniform subsonic crossflow from a flat 
plate. The jet-to-crossflow velocity ratio must fall between three and 
ten, A flow chart outlining the structure of the JASIB76 module is 
shown in Figure 6. 

The JASIB76 executive program oversees the computation procedure by 
systematical 1 y cal 1 i ng each of the subrouti ne groups • Contai ned wi thi n 
the executive program, in the form of comment cards, is a brief history 
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Figure 6. - JASIB76 Flow Chart. 
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of the JASI programs and a description of each subroutine within the 
module* The executive program In the module JAS 1676 receives the 
FORTRAN CAU Statement from the outside program, In this case WBABL* 

The CALL statement to JASI676 contains the arrays specifying the field 
points of Interest and the jet perturbation velocity arrays associated 
with the field points* Both the field points and velocity components 
are In the panel code reference coordinate system* Variables passed to 
JASIB76, from the panel code, through COMMON statements are the 
following: the location of the Jet center; Jet diameter, 0; Jet to 
crossflow velocity ratio, R; Jet injection angle, 6; angle of attack, a; 
yaw angle, g; and the Jet roll angle, y* The dimensional units given in 
the list above are all in the panel code reference coordinate system* 

The subroutine sections called by the executive program are described 
below in the order of their calling, 

DATA. This subroutine reads from the input file the options and 
parameters for the jet/ flat plate model. The values In the parameter 
sets are dependent on the following jet properties: velocity ratio, R; 
Injection angle, 5; jet exit shape and profile, 

SEGCL, Thi s section approximates the jet plume center! ine 
geometry. The Jet centerline Is broken Into linear, constant strength 
source segments to model the Jet entrainment. The SEGCL subroutine 
calls six other subroutines in order to divide the jet centerline and 
determine the source strength on each segment. 

SEGVC. The SEGVC section approximates the geometry for the vortex 
paths. As In the SEGCL section the vortex paths are segmented and 
constant strength vortices are placed on each segment. Only the 
starboard vortex path Is approximated; the geometry of the vortex paths 
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Is synifnetrleal about the x-z plane In the Jet coordinate system* The 
vorticity on the port vortex segments Is equal In strength and opposite 
In sign to that of the starboard segments. 

TOTVVC, This section calculates the velocity Induced at the field 
points by the contrarotating vortex system. Contained in TOTVVC Is the 
calculation of the transformation matrix for relating the jet coordinate 
system, as shown In Figure 2, to the panel code reference coordinate 
system, as shown In Figure 1. The transformation matrix Is used to 
transform the field points to the Jet coordinate system and the velocity 
components from the Jet coordinate system Into the panel code reference 
coordinate system. The subroutine VVC, which Is called by TOTVVC, is 
the routine which actually calculates the perturbation velocities for 
each field point. 

TOTVCL . The TOTVCU subroutine, which uses subroutine VCL, Is 
synonymous in function to TOTVVC, In that It calculates the perturbation 
velocities at the field points due to the jet plume entrainment. 
Subroutine VCU Is called by TOTVCL. VCL calculates the perturbation 
velocity for each field point. 

WAKSlIB. Subroutine WAKSUB is used to model the wake region 
downstream of the jet In a crossflow. At the time of this writing the 
wake subroutine Is Inactive in the JASIB76 program module. The reason 
for the lack of a wake subroutine Is that the only wake model available, 
at the time the program was coded, was an empirical correction 
subroutine which works directly with the surface pressures; only the jet 
perturbation velocities are availabVe in the jet module. At present 
wake modeling is handled by a post-processor unit, which modifies the 
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pressure distribution downstream of tlie jet after the tnain program has 
finished execution* 

Output from Che JASI676 program Is provided through three 
channels: the CALL statement addressing JASIB76. COMMON blocks, and the 
panel code output file* The jet perturbation velocity components for 
each field point are provided through the CAU statanent; the 
perturbation velocities are added to the values provided by the CALL 
statement at the Initiation of JASIB76. The Isentropic jet thrust 
components, In the panel code reference coordinate system, are output 
through the COMMON blocks which link the panel code to the Jet module* 
The output, written to the output file, IS dependent on the diagnostic 
and output options specified in the jet Input parameters. Diagnostic 
output from 0ASIB76 consists of the following; all Input parameters and 
control point locations; geometry for the jet centerline and the 
starboard vortex; vortex induced velocities at a field point for all 
vortex segments; total vortex induced velocity for all field points; 
entrainment Induced velocities at a field point for all centerline 
segments; total entrainment Induced velocity for all field points; and 
the jet thrust components* 

Recall from the previous description of WBABL that the overlay 
WBAERO occupies the largest quantity of core memory, and that the 
JASIB76 code must be placed in this overlay. For this reason the core 
memory requirements of JASIB76 are reduced. The primary method of 
reducing the core memory requirements of a program executing on the 
COC-7600 computer 1s to place as many variables as possible In the Large 
central Memory (LCM) region of the computer. All variables In labeled 
COMMON arrays within JAS 1876 are placed In LCM, As a result of placing 
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dll possible arrays in LCM, the JASI676 program requires only three 
variable arrays to be in the core memory. These are the perturbation 
velocity components. Each of the velocity arrays Is dimensioned to the 
number of control points* The only other arrays In core memory of any 
significant size, are the field point locations, and they are required 
to be In core memory by the panel code, WBABL. 

The Combined Gomouter Code. VIBVIJAS 
The combined computer code, known by the acronym WBWJAS, is formed 
by adding the JASI876 module to the MBABL panel code as a set of 
subroutines In the MBAERO overlay. A flow chart of the WBWJAS execution 
is shown in Figure 7. The flow chart shows that the structure and flow 
of the WBWJAS program is identical to WBABL, except for the addition of 
Jet effects before the potential flow calculation. If the WBWJAS 
program is executed without jets in crossflow, the execution and output 
would be identical to that of the WBABL program. 

The execution of WBWJAS with a jet in crossflow will differ from 
WBABL as follows. After lofting the geometry for the body and wing 
panel networks, the program reads to determine If there will be any jets 
in crossflow and if so how many. The panel code next determines the 
velocity field perturbations due to unit singularity strengths on the 
source and vortex panels. This is the same as WBABL. When all geometry 
is known, WBABL calls the JAS1B7$ module, passing It the field points 
and jet perturbation velocities for each field point* When execution 
returns from JASIB76, the value of the jet perturbation velocities is 
the input value plus the value determined for the jet that was just 
addressed. After all jets in crossflow have been addressed, the panel 
code treats the jet perturbation velocities as if they are a part of the 
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freestream flow. The right side of the matrix equation, the boundary 
eondltionst Is computed using the jet-1 n-cross flow perturbation 

f 

velocities as components of the freestream flow. The matrix equation Is 
solvedt as It would be In the non-jet in crossflow casei for the 
singularity strengths* With the slngularivy strengths known* the 
velocity at each control point is the sum of the freestream velocity and 
all perturbation velocities* The pressure coefficients at each control 
point are then found using Bernoulirs equation* This completes the 
potential flow calculation for WBWJAS* 

The streamline tracing and boundary layer calculation sections 
follow the potential flow calculations* These sections do not have to 
be exercised when executing the panel code and for this thesis they are 
not exercised* The final computation In the panel code Is the numerical 
integration of the pressure distribution in order to determine the 
forces and moments* The decision block In the flow chart (Figure 7) 
following the force and moment calculation, ITR»ITRMAX, refers to the 
number of angle of attack cases being Investigated* If only one case is 
used the program would stop execution after the first complete run. 

The WBWJAS program is not without limitations. The program 
contains most of the limitations of the WBABL and JASIB76 programs plus 
others* The reader Is reminded of three main limitations to JASIB76; 

It models a round uni form, subsonic jet Injected perpendicularly Into a 
uniform subsonic crossflow; the jet-to-crossflow velocity range Is valid 
from three to ten; and the model properties at this time are based On a 
jet/flat plate configuration* The limits of WBABL are as follows; It 
does not model the interference effects of jets in crossflow; it is 
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Hmiteii to 1500 total control points on one side of the X-l symmetry 
plane; and It is limited to subsonic attached flow only. 

The limitations of the combined program, W5WJAS, which are not 
mentioned above are included now. The code is restricted to 
approximately 660 total control points on one side of the X-Z symmetry 
plane, due to memory restrictions. The influence of the lifting 
surfaces on the jet properties are not accounted for. There are no 
interactions between multiple jets in crossflow, and the wake region 
downstream of the jet core is not modeled. 

The input/output structure of WBWJAS 1$ simply a combination of the 
U6ASL and JASI876 codes, with all options of each code still available. 



SECTION IV 

CONFIGURATION TESTING 

Three basic conflguratlgns are tested using the WBWJAS paneVcode ■ 
with Jet-1 n«cross flow effects* The three configurations are a finite 
flat plate, a body of revolution, and a symmetrical airfoil section. 

The panel code configurations are paneled to match as closely as 
possible the configurations which were tested in low speed wind 
tunnels. The first configuration tested, the flat plate, Is used as a 
tool to verify the WBWJAS program execution. The body of revolution and 
the symmetrical airfoil are both examples of easily modeled generic 
configurations, for which wind tunnel pressure data is available. 

Flat Plates 

Three flat plate configurations are utilized in the development of 
the WBWJAS program. The first is used only for program development work 
to verify any modifications. It is a thin, rectangular grid, eighty- 
four control point flat plate with no thickness. The second flat plate 
is also used for developmental testing and is not presented as a data 
source. This plate is rectangular with a finite thickness. The upper 
surface is a radial grid with the jet located at the center. The lower 
surface Is a rectangular grid; Its edges are level with the upper 
surface. The final flat plate, which is used for actual testing 
purposes, is a thin, rectangular grid plate. The semi -span dimension Is 
iOO units with the chord length variable. The fi nal chord 1 ength is 165 
units and the panels are five units square. These dimensions yield a 
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final configuration of 660 panels* The jet In crossflow Is located at 

t 

(75t0i0) and has a variable diameter* The final jet diameter Is ten 
units* As tested the flat plate dimensions In Jet diameters are length 
16*S 0, width 20 D; the jet Is located on the centerline 7*5 0 from the 
leading edge* The reference coordinate system Is located at the leading 
edge along the plane of symmetry* The reference angles a, 8* and y are 
all zero* indicating the jet coordinate system Is parallel to the 
reference coordinate system* 

The flat pi ate. configuration 1$ used to determine the maximum 
number of panels a configuration can have* and to determine the 
acceptable range of panel densities* which give satisfactory results* In 
the near jet region. The maximum nunber of panels Is found by varying 
the plate chord. The acceptable range of panel densities can be found 
by varying the jet diameter. 

The wake region downstream of the jet core Is modeled by an 
empirical wake correction scheme as described In reference 8. The wake 
correction Is In a post-processor unit to the panel code and modifies 
only the surface pressure distribution. The wake cc.^rectton uses the 
empirical pressure distribution, along the downstream ray from the jet 
orifice, to estimate a correction factor which is applied to the 
pressure distribution within the wake region. 

Body of Revolution 

Hany^ experimental Investigations of jet-in-crossflow flow 
properties have placed attention on flat plate configurations, or 
jet-ln-li ftlng surface configurations. These configurations have little 
1 f any curvature of the 11 fti ng surface near the jet 1 n crossfl ow . Few 
Investigations have been made of a jet/aerodynamic-surface configuration 
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with a high dagree of surface curvature near a jet tn a cross flow (e*g«. 
a jet/fuselage conMguratlon), A single engine v/STQL aircraft will 
probably have lift jets located In the fuselage, for the simple reason 
that locating them elsewhere would require complex ducting, adding extra 
weight and expense. The simple generic model for a fuselage Is a body 
of revolution. A recent wind tunnel Investigation.^^ at NASA Ames 
Research Center, studied the surface pressure distribution on a body of 
revolution with various jet configurations. The body-of-revolutlon 
configuration studied at NASA Is modeled here with the U6WJAS computer 
code, for the case of a single jet In a crossflow with perpendicular jet 
Injection. 

The body-of-revoluti on model tested by NASA Is shown In Figure 8, 
which is reprinted from the test repprt.^^ The model is four inches in 
diameter and eighty-one inches in length. The jet diameter Is 
1.94 inches. The model has a wooden, streamlined nose, and Is strut 
supported at the rear. The jet center for the single jet case Is 
located twenty-two inches from the nose. Two jet nozzle exits were 
used. The first was flush with the surface of the body, which is 
referred to as a contoured exit. The second, a flat top exit, was level 
with the top of the model. 

The computer configuration of the body of revolution Is paneled to 
approximate the control point locations to the pressure ports on the 
wind tunnel model* The computer configuration is divided into six axial 
sections known as patches. The patches are paneled using the 
body-of-revoluti on geometry option in WBWJAS. The patch structure 1s 
shown tn Figure 9, along with an end view showing the panel strxture in 
the jet region. The nose, forebody, aftbodyi and tall patches are 
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coarsely paneled, such that they maintain the body shape with a 
relatively small number of panels. The jet patch Is simulated with a 
fine panel grid, approximating the pressure port locations on the test 
model. The jet patch contains 540 panels and all other patches combined 
have 120 panels. The reference coordinate system origin for the body- 
of-revolutlon model Is located ’at the nose of the model, as shown In 
Figure 9. All length units are relative length units based on a body 
diameter of four units. The jet In crossflow Is located at (22,0,2.0) 
for the flat top exit, and (22,0,1.732) or (22,0,1.5) for the contoured 
exit. The jet Is tv (0 units In diameter. The orientation angles a, 
g, and Y are all zero. 

A wake correction is not applied to WBWJAS for the body-of- 
revoluti on models. As will be seen In the next section, the data shows 
that a wake correction Is not needed for this configuration. 

Symmetrical Wing Section 

Throughout the early sections of this thesis it is stressed that 
WBWOAS does not model the Influence of lifting surfaces on the jet-in- 
crossflow properties. To investigate the li fting-surface influences, a 
symmetrical wing section is modeled. The symmetrical wing section 
modeled is similar to a jet/wing configuration that was tested by 
Mlkolowsky,^® in the Georgia Institute of Technology nine-foot low-speed 
wind tunnel. The data from the W8WJAS model are compared with data from 
the wind tunnel model. 

The wind tunnel model spanned the entire test section with a 
NAC A 0021 airfoil section whose chord was 15.37 inches. Two jet 
configurations were tested and are now described. The first jet 
configuration had a diameter of 1,5 Inches, with the jet centerline exit 
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location at either 25, 45, or 65 percent chord. The second had a jet 
exit diameter of 3*0 Inches exiting at 45 percent chord. The jet-to- 
crossflow velocity ratio was varied In the range from two to twelve. 

The tests were performed at three wing angles of attach: 0®, 6®, and 
9®. 

The computer configuration paneled for UBUJAS Is that of a finite 
semi-span wing. It has a NACA 0021 airfoil section with a chord length 
of 15 units, and an aspect ratio of twenty. The jet diameters used are 
1,5 units and 3,0 units, and are placed In the wing at identical 
locations to those on the wind tunnel model. The configuration is 
tested at jet-to-crossflow velocity ratios of four and eight, and angles 
of attack of 0® and 6®. 

The wing is not paneled such that the control points coincide with 
the pressure port locations. The wing model Is divided Into three 
patches. They are a near-fteld patch, a mid-field patch, and a far- 
field patch. The near-field patch extends from the symmetry plane, on 
which the jet Is located, to a point 4.5 units In the spanwise 
direction. The near- field patch contains ten spanwise sections, each 
containing forty-six panels (thirty on the lower surface and sixteen on 
the Upper surface). The mid- field patch, extending from 4*5 units to 
15.0 units In the spanwise direction, contains four spanwise sections. 
Each spanwise section contains sixteen panels on the lower and upper 
surfaces, for a total of thirty-two per section. The far- field patch 
extends from 15,0 units to 150,0 units In the spanwise direction, and 
contains a total of forty-eight panels In six spanwise sections. 

The jets are always located along: the plane of symmetry and exiting 
the lower surface of the wing perpendicular to the freest ream flow. In 
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the wtnd tunnel model the jets were perpendtcul ar to the chord line, 
giving non-perpendicular jet Injection angles. For the computer model 
the Injection angle is always perpendicular to the crossflow, regardless 
of the wing angle of attack. The jet is located on the wing lower 
surface. The reference angles are the following: a, equal to wing 
angle of attack; s> always zero; y» always one hundred eighty degrees, 
so the jet exits from the wing lower surface,- 

A wake correction is made to the computer generated data. This 
wake correction is Identical to the one used on the flat plate 
configuration. This wake correction is applied only to the surface 
pressure data, and not to the force and moment data. As a post 
processor section, the wake correction scheme could not be supplied with 
the necessary data to correct the force and moment data. 



SECTION V 

RESULTS AND CONCLUSIONS 


Flat Plates 

Flat plate models are used primarily to debug the program and 
validate the numerical methods* The third fiat platei mentioned earlier 
as having a variable plate chord and variable jet diameter, is used to 
find the maximum number of control points permissible. The maximum 
number of control points that can be used In the N8WJAS program Is 
dependent on the amount of core memory available, and is independent of 
the configuration being tested. The maximum number of control points Is 
found by varying the flat plate chord length while maintaining a 
constant panel size. For the CDC-7600 computer system at Ames Research 
Center the maxlmiim number of control points Is six hundred sixty. 

The jet diameter Is varied on the flat plate, while maintaining a 
constant panel size, to determine an acceptable range and value, for the 
ratio of the jet diameter to the control point spacing, Rg/sp^ The 
ratio Rgy^p Is synonymous to the panel density, and is Important only 
within the near-jet region, i.e., within approximately four jet 
diameters of the jet orifice. The criterion for finding the acceptable 
range of the ratio Rg/sp is the agreement of surface pressure contour 
co-plots. This criterion can be interpreted as finding the minimum 
number of data points per unit area necessary for the particular contour 
plot program being used. This is not the Intent of the criterion. The 
intent of the criterion Is to find the acceptable range for Rg/gp 

42 



43 


required to yield an acceptabie pressure distribution on the surface of 
interest* A fnaxlmutn ratio can be found by testing for the effects of 
roundoff error in either the output data or resulting contour plots* 

Good agreement between the computer pressure contours and the 
experimental pressure contours for a jet/ fiat plate model are found 
for 1*0 < Rq/sp < 5*0 * At values of Rp^gp greater than 5*5 the effects 
of roundoff error become evident, and below 1*0 the data is too sparse 
for the contour plotting program to produce smooth contour Tines. 
Agreement of pressure contours is found to begin at values of R^^gp = 

2*0 and continue to 4.0 • Thus the value of Rg/sp used for the flat 
plate models is %/sp * 2.0 • 

Surface Pressure Data 

A comparison of pressure coefficient contour Vines is made, for 
three di^-ferent sets of data, in Figure 10. The three data sources are 
all for a uniform round jet Issuing perpendicularly into a uniform 
crossflow from a flat plate. The three data sources are experimental 
data from reference 25, WBMJAS, and JASIB^ which models the flat plate 
as an Infinite plate by the method of Images, The figure shows 
satisfactory agreement between the three data sources for all contour 
lines shown. 

Force and Moment Data 

The interference lift and pitching moment effects of the jet in a 
crossflow on the flat plate are also documented. The interference lift 
is npn-dimenslqnalized by the isentropic jet thrust and is defined 

^ ■■ ^jet on * *~jet off ~ ^j 

' Y ' 
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The Interference pitching moment about the Jet center Is non- 
dimensional 1 zed by the product of the Isentroplc jet thrust and the jet 
diameter, and Is defined 

^.1 ",i .1»t on ~ let off 

V V 

The lift and jet thrust for the flat plate are defined to be positive In 
the negative "Z" direction of the panel code coordinate system* The ■ 
pitching moment is defined positive In the negative "Y" direction for a 
right hand coordinate system. 

The Interference lift and pitching moment on the flat plate are 
given in Table 1* The data from the MBWJAS program Is presented with 
and without a wake correction. The data from reference 8 for the JASIB 
program, with and without a wake correction, Is provided for comparison 
with MBWJAS data. It should be noted that the forces and moments 
provided by JASI8 and the experimental data are produced by numerically 
Integrating the surface pressure distribution over a reference circle. 
This reference circle Is centered at the jet center and has a diameter 
of eleven jet diameters. Force and moment data from UBWJAS are 
determined by using an area which approximates the reference circle for 
numerical integration of the surface pressures. The differences In the 
force and moment data reflect the differences In the surface pressure 
contour plot, Figure 10, for this flat plate configuration. 



TABLE 1. - INTEKFERENCE FORCES AND HOHEHTS 


CASE 

HODEL 

’i V 

EXPERIHENTAl 

Tj Tjd 

FLAT PLATE* 

MBWJAS Data, no wake 
Data, no wake 

HBViJAS Data, >»ke correction 
JASIB Data , wake correction 

•0.20 -0*091 
-0.20 -0.07S 
-0,23 0.012 
-0.25 0.018 

^0.25 0.013 
-0.25 0.013 
-0.25 0.013 
-0.25 0.013 

BODY OF REVOLUTIOI^ 

Flat Top Exit, R=8, ZJCT=2*00 
Contoured Exit, R=8, ZJCT=i.732 
Contoured Exit, R=B, ZJCT»1.50 
Contoured Exit, R=4*7, ZJCT-1.Z32 

-0.04 0.12 
-0*10 0.09 
-0*16 0.06 
-0.26 0.17 



^Expeirimental data fram Reference 25. 

^No experimental force and moment data available* 
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Body of Revolution 

The t}Ody-of-revo1ut1on model Is Intended to Investigate the 
usefulness of the UBUJAS program for modeling bodies with high degrees 
of curvature near the jet In a crossflow*. The model utilized, as 
described earlier, has a body>to-jet d1air«ter ratio of two. The body* 
of*revolut1on data 1$ presented without a wake correction. Data Is 
presented for four body-of-revolutlon configurations. The first three 
configurations have a common value of the jet-to-crossflow velocity 
ratio of eight; the jet center location “.Z" value Is varied In the three 
configurations. The fourth configuration Is with the velocity ratio 
equal to 4.7 • 

Surface Pressure Data 

The surface pressure data, provided by WBWJAS, are compared with 
the results of the NASA experlment.^^ The data is presented as 
interference pressure coefficients. These are defined as the pressure 
coefficient with the jet on, less the pressure coefficient with the jet 
off. Both values are measured at equivalent test conditions. The data 
Is presented as plots of axial variations In interference pressure 
coefficients, aCp. The axial plots are presented for constant Y/0 cuts 
of 0,0, 0.24, 0.47, and 0.84. 

Figure 11 presents four surface pressure plots for a jet center "V 
value (ZJCT) of 2.0 • This configuration Is representative of a flat 
top jet exit. Experimental data for a similar configuration is also 
shown In these plots. The data from WBWJAS, plotted as a solid Tine, 
agrees well with the experimental data. 

Figures 12 and 13 present data for ZJCT values of 1.732 and 1.50, 
respectively. These values of ZJCT are Intended to model the contoured 
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jet exit* The value of ZJCT > 1*732 Is chosen because the experimental 
model has jet side walls to a height of 1,732 Inches above the body 
centerline* The ZJCT ■ 1*50 value Is chosen to Investigate any trends 
that may develop by submerging the jet within the body. The pressure 
plots In Figures 12 and 13 show that the WBHJAS data matches closely the 
experimental data for a contoured jet exit bn all but the symmetry plane 
axial cut* The UBWJAS program does not predict the negative pressure 
field to the aft of the jet along the Y/D ■ 0*0 axial cut. The program 
shows a deceleration of the flow forward of the jet orifice and an 
acceleration region away from the jet orifice, downstream of the jet. 

The experimental data, on the other hand, shows no deceleration of the 
flow forward of the jet, and a definite low pressure region aft of the 
jet, i,e,, a small wake region. The data of Figure 13, for ZJCT • 1,50, 
Indicates that as the jet is submerged Into the body, the computer code 
does not predict the low pressure region aft of the jet, along the Y/D « 
0*0 cut. The submerging of the jet-1 n-crossflow exit does slightly 
Improve the data for off centerline cuts (Y/D * 0,0), when modeling the 
contoured exit jet* 

Figure 14, which presents data for the contoured jet exit with the 
velocity ratio equal to 4*7, shows similar trends to the R » 8 contoured 
jet exit configuration. This Indicates that discrepancies between the 
experimental and computer generated data are Independent of jet-to- 
crossflow velocity ratio, 

A wake correction Is not used bn the body-of-revolutlon models. 
Close agreement between the experiment and computer generated data for 
all axial cuts except the contoured exit symmetry plane axial cut 
Indicates a wake correction Is not necessary. This author believes the 
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absence of an extensive wake region aft of the jet In a crossflow Is a 
result of surface flow around the body Into the region aft of the jet 
orifice* This surface flow Is possible due to the high degree of body 
curvature which results In little actual body surface area within the 
region aft of the jet orifice* 

Force and Moment Data 

The Interference lift and pitching moment data for the body-of- 
revolution models Is presented In Table 1* The interference lift and 
pitching moment are defined the same as for the flat plate 
configurations. Wind tunnel experimental data 1$ not available for the 
body-of-revol utl on configurations* 

The Interference lift data shows a lift loss due to interference 
effects, which is to be expected* The steady increase of the lift loss 
Is also expected from the submerging of the jet Into the body, or the 
reduction of the jet*to-crossflow velocity ratio. The nose-up pitching 
moment is expected, although the trend of decreasing nose-up pitching 
moment with jet submersion Into the body Is not expected. Based on 
pitching moment changes for flat plate configurations, one would expect 
that as the Interference lift loss Increases, so would the Interference 
pitching moment, although for this case, the majority of the lift loss 
is not aft of the jet. The Increase In Interference pitching moment Is 
expected for the fourth case (8R 4), relative to the second (BR 2), 

Symmetrical Wing Section 

A symmetrical wing section modeV is used to Investigate the 
accuracy of the jet/ flat plate model where simulating a jet/wing 
configuration. The symmetrical wing section model used, as described 
earlier, 1$ a NACA 0021 wing section with a fifteen unit chord and an 
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aspect ratio of twenty. The surface pressure data Is presented with a 
wake correction which 1$ similar to the flat plate wake correction. Six 
configurations are tested; a base configuration and five others which 
are each a single parameter perturbation of the base configuration. 

The wing section modeled by the computer program Is a three- 
dimensional semi-span wing, whereas the experimental model Is a two- 
dlmenslonaT full-span wing. In order to compare the two sets of data, 
it Is necessary to adjust the computer model’s span such that three- 
dimensional (tip) effects within the region of interest are 
negligible. The tip effects are found to be Insignificant, within the 
region of Interest, by adjusting the full -span aspect ratio to twenty. 
This places the wing tip one hundred jet diameters from the jet orifice; 
the region of interest for interference pressure plots extends to eleven 
jet diameters from the jet orifice. 

Six different wing configurations are modeled with WBWdAS. The 
wing configuration parameters which are varied are the following: the 
jet exit location In percent chord, Xj/c; jet diameter in percent chord, 
Dj/c; jet to crossflow velocity ratio, R; and the wing angle of attack, 
cu The base configuration has the following parameter values; Xj/c » 
45%; Pj/c = 10%; R » 8; a* 0.0®. Changes are made to the base 
configuration as follows; Xj/c to 25% and 65%; Oj/c to 20S; R to 4; and 
a to 6®. The data obtained from each computer configuration Is compared 
with experimental data and presented as surface pressure plots and force 
and moment data. 
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Surface Pressure Data 

Data from each configuration is presented in the form of an 
Interference surface pressure (aCp) contour plot. Each graph Isa 
co-plot of data from the UBUJA5 computer program and experimental 
contour from wind tunnel tests.^® The plots show the Interference 
surface pressure coefficient contours on the lower surface of the 
wing. The pressure plots are given In Figures 15 through 20 with the 
base configuration shown In Figure 15. Each change to the base 
configuration Involves only one parameter each and are shown as 
folVws; Figure 16 Is for a = 6® ; Figure 17 for the forward jet 
location* Xj/c >25%} Figure 18 the aft jet location, Xj/c » 65%; 

Figure 19 is the lower value of jet-to-crossflow velocity ratio, R » 4; 
and Figure 20 Is the large diameter jet, Oj/c = 20%. Each aCp contour 
line is labeled If space permits; If a line is hot labeled Its value Is 
that of the closest neighboring contour line. 

The following observations can be made about all of the contour 
plots: in the near jet region, less than two and a half jet diameters 
from the jet orifice, the interference surface pressures are predicted 
satisfactorily by the WBHJAS program; for the remainder of the region of 
Interest the WBWJAS program is predicting stronger jet interference 
effects than those observed in the experimental data. The differences 
between the two sets of data might be attributed to ^t^^ lack of lifting 
surface Influence on the jet plume properties^ 

The effects of changing angle of attack to ot = 6® do not appear to 
be strong. One should remember that the jet is on the wing lower 
surface. The effects may be more pronounced if the angle of attack had 
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been changed to a* -6^« Experimental data Is only. available for 
positive angles of attack. 

The effects of jet movement fore and aft on the wing can be seen In 
Figures 17 and 18. Moving the jet forward on the wing tends to magnify 
the disagreement In the two sets of data. Moving the jet aft on the 
wing tends to decrease any disagreement In the two sets of data. These 
two observations may be due to the greater lifting surface effects near 
the leading edge of the wing. 

The effects of decreasing jet-to-crossflow velocity ratio (R) are 
seen In Figure 19. This decrease In velocity ratio Is similar to 
stronger lifting surface Influences. Again the lack of lifting surface 
Influence on the jet In crossflow may cause the computer model to 
predict stronger jet-1 n-crossflow Interference effects. 

The effect of Increasing jet diameter, as shown In Figure 20, 
appears only to magnify the results for the base configuration. The 
jagged appearance of the experimental data contour lines may be due to 
sparse data. It has been this author's experience that sparse data 
often lead to Inaccurate contour lines. 

Force and Moment Data 

Interference lift and pitching moment data are not presented here 
for the symmetrical wing. The reader Is reminded that the force and 
moment data from VjBWJAS is only available without a wake correction, 
therefore, It Is not representative of the real flow field. The author 
of reference 26 Indicated that he encountered balance problems during 
the wind tunnel tests, thus the experimental data is not reproduced 
here. 
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General Conclusions 

Results from the configurations tested show the jet/aerodynamlc- 
surface model* as presented here* Is capable of approximating the 
adverse aerodynamic interference effects associated with some generic 
aerodynamic configurations in transitional flight* The results of the 
flat plate models indicate that the jet/aerodynamic-surface model 
presented here Is correctly approximating the surface pressure 
distribution on a flat plate model when a wake correction is used in the 
viscous wake region* The body-of-revolution model demonstrates that a 
wake correction Is not always necessary* and that the Jet interference 
effects on this body are similar to those on the flat plate model* The 
surface pressure data from the symmetrical airfoil models indicate the 
lifting surface effects on the jet In crossflow are not detrimental to 
this Jet/aerodynam1c>5Urface m'^dei* These favorable results should 
provide the necessary justification to continue development on the 
iterative method proposed* 

Continued Research 

As a close to this report the author would like to suggest the 
following areas for research and development on the WBWJAS program: 
development of a watce model which manipulates the jet perturbatibn 
velocities In the wake region; extension of the Jet-in-crossflow model 
to handle varying jet Injection angles; extensions to handle various jet 
exit velocity profiles and shapes; Increases in the range of valid jet- 
to-crossflow velocity ratios; and most Important, to make' the program 
iterative for the effects of lifting surfaces. 
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